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Abstract 

We study the role of diamagnetic effects on the transport properties of metal- 
lic magnetic multilayers to elucidate whether they can explain the Giant 
Magnetoresistance (GMR) effect observed in those systems. Realistic Fermi 
surface topologies in layered ferromagnets are taken into account, with the 
possibilities of different types of orbits depending on the electron spin. Both 
configurations, with ferromagnetic and anti-ferromagnetic couplings between 
magnetic layers, are considered and the transmission coefficient for scattering 
at the interface boundary is modelled to include magnetic and roughness con- 
tributions. We assume that scattering processes conserve the electron spin, 
due to large spin diffusion lengths in multilayer samples. Scattering from the 
spacer mixes different orbit topologies in a way similar to magnetic break- 
down phenomena. For antiferromagnetic coupling, majority and minority 
spins are interchanged from one magnetic layer to the next. Cyclotron orbits 
are also traveled in opposite directions, producing a compensation-like effect 
that yields a huge GMR, particularly for closed orbits. For open orbits, one 
may get the inverse magnetoresistance effect along particular directions. 

PACS numbers: 72.15.Gd, 75.70.Pa, 73.50.Jt 
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The Giant Magnetoresistance (GMR) effect discovered in magnetic multilayerslll is known 
to have many potential applications, including recording materials! and new promising elec- 
tronic devices based on spin effects.!!! This field of magnetism, that now involves surface 
physics, metallic films, magnetic small particles and nanostructures, is currently very active 
in both, basic and applied research. A typical multilayer system consists of ferromagnetic 
films (usually iron or cobalt), separated by thin non- magnetic metal layers. The GMR ef- 
fect appears when the ferromagnetic films are coupled through the spacer in an antiparallel 
configuration. The application of an external magnetic field yields a negative magnetore- 
sistance, which is typically of the order of 50~60% for the Current In-Plane (CIP) set 
up. The effect is attributed to the existence of spin dependent electronic scattering!!, in 
such a way, that transport properties can be described in terms of a two current model, 
with two spin channels^. Further developments include modelling of the interface roughness 
and quantum effects in very clean samples!. Within the semi-classical approach, most of 
the calculations in the literature integrate the Boltzmann equation with different boundary 
conditions, following the original proposal by Fuchs and Sondheimer for thin filmsiH. Dia- 
magnetic effects of the electron orbits have been neglected when compared to the interface 
scattering, which is consider as dominant in most contributions. We would like to argue in 
this paper that this is not a sensible approximation. In a ferromagnetic metal, like iron or 
cobalt, there is a huge internal magnetic field, resulting in very different Fermi surfaces for 
each spin§h0. As a consequence of the large spin splitting, we may have different electron 
occupations at the Fermi level for each spin band, and depending on the direction of the 
internal field, we may vary the orbit topology in a different manner for each spin. This is 
the case of iron, where the minority Fermi sheet supports an open orbit along the < 001 > 
direction. In contrast, the majority Fermi surface displays a multiply-connected open orbit 
in the < 110 > direction!. A recent band calculation for ferromagnetic fee cobaltS shows 
a variety of situations: both spins present closed orbits for the (001) plane; on the other 
hand, for the (210) plane, the majority-spin orbits are all closed while the minority Fermi 
surface supports open orbits. Most of the scattering mechanisms at the spacer, which is 
paramagnetic or non-magnetic, conserve the spin (no spin-flip). With the present state of 
the art in fabricating multilayers, the spin diffusion length for spin polarized carriers are 
of the order of the sample dimensions!. This means that, whether transmitted or reflected 
back at the interfaces between layers, the carriers remains in the same spin Fermi sheet. For 
the anti-parallel configuration, the internal magnetic field is reversed when going from one 
magnetic layer to the next, with the corresponding interchange between the majority and 
minority spin Fermi surfaces. Similar effects has been considered by Cabrera and FalicovEi, 
to study transport near magnetic domain walls. The same concepts lead to the following 
scenario: i) a different availability of states for each spin at the Fermi level, producing a 
spin valve effect when the electronic current flows through the multilayer structure; ii) due 
to the strong internal field, diamagnetic effects are important. Since the Fermi surfaces 
are different for the two spin states, scattering from the spacer couples electron orbits of 
different forms and/or topologies at both sides of the interface, in a way similar to magnetic 
breakdown phenomena^. As usual, diamagnetic effects are important when x = 3> 1, 
with uj c being the cyclotron frequency and r the bulk mean collision time. 

In this paper, we propose a simple model that incorporates the above affects, with the 
presence of the sample walls as an essential ingredient that may affect the connectivity 
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of the cyclotronic orbits. If the mean collision time for scattering from the interface is 
much smaller than the cyclotron period T c , electrons at the Fermi surface have no time to 
explore the orbit topology, even if the condition % ^> 1 is satisfied. The interesting case 
corresponds to the situation when both time scales are of the same order, which implies mixed 
orbit topologies. We note that at high magnetic fields, the qualitative magnetoresistance 
behavior does not depend on details of the Fermi surface, but only on its topologyB. Simple 
spherical shapes are then used in our numerical calculation to model closed and open orbits. 
The sample geometry adopted in our model, is that of a three-layer sandwich, with wide 
ferromagnetic layers separated by a thin non-magnetic spacer. In a simplified picture, the 
spacer is considered as an interface between magnetic layers, whose roughness is modelled 
through a potential barrier Vq, independent of the spin. In addition, due to the exchange 
splitting Vi j the electrons are subjected to different magnetic potentials, depending on their 
spin. We have calculated transport properties for the two magnetic configurations, the one 
with anti-parallel coupling between magnetic layers (APC), and the other one with parallel 
alignment (PC), and measured their difference. We assume that the sandwich is repeated 
periodically, to generate the multilayer structure. 

To calculate the conductivity tensor, we use the Chambers' path integral approach, which 
is particularly suited to treat the contribution of entangled networks of electron orbitsS. 
Within this scheme, the conductivity is given by 
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Stc^Jfs |v A (k)| dS ' (1) 
with v A (k) being the velocity vector for an electron on the Fermi surface, and 

AkA(r) = T (k) v^exK * ~ * 0(k V , (2) 

J— oo T 

being the effective mean free path for wave vector k and spin A. The above formulae have 
been established using the Chambers kinetic approach^, which is equivalent to solving the 
Boltzmann equation, taking into account terms to all orders in the magnetic field. As usual, 
the time dependent v kA (t) = (1/h) V^e A (k) in (g), represents the velocity of an electron 
which is in the state k at time to(k), and is obtained solving the equation: 



hk = -— v k x B . (3) 



The integration of (|3|) yields a tt/2 rotation between the orbits in real and fc-space, apart from 
a scale factor that depends on the magnetic field. We choose the perpendicular 
to the interfaces, with the internal magnetic field being B = Bz or B = —Bz, depending 
on the kind of magnetic coupling between ferromagnetic layers. The cyclotron orbits are 
then realized in the xy-plane, and we are just calculating components of the relevant tensors 
projected on that plane. We take a complete network of orbits and divide it into n suitable 
pieces named a, (3, to describe the proper Fermi surface topology. Further, we label the 
state k of an electron by a piece-number and a time t. Electrons moving on a particular piece 
orbit a are in the states (t, a), with < t < t a . We assume that at some points on each orbit 
piece, due to the interface scattering, there is a probability for an electron to be transmitted 
to a different piece orbit or to be reflected back, remaining in the same Fermi surface sheet. 
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For simplicity, we have considered the same time t a for all the pieces (isochronous pieces), 
thus t a is the elapsed time for electrons, moving on the a-piece, between two scattering 
processes from the interface. Different networks, with different connectivities, are generated 
depending on the time scales of the problem, i.e. depending on the parameter 77 = t a /T c , 
where T c is the cyclotron period. 

With this notation, we may write relation (pD in terms of a discrete set of quantities 
|A( a )|, in the form 

A(t, a) = A (a) e- t/r + e~ t/r [ v(t', a)e t,/r dt' (4) 

Jo 

where a = l,...,n. The A( a ) terms take all the information about the past history of 
electrons on the network before reaching the a— piece at time t = 0. In the limit case of 
total reflection from the interface, electrons remain in pre-defined orbits for each spin state. 
On the other hand, in the limit of total transmission, the electron moves in a complex, but 
coherent, network of orbits mixing the Fermi topology of the magnetic layers. So, one gets 
a kind of artificial 'breakdown' interface-driven mechanism between this two limits. Further, 
we include all the information about the connectivity among the pieces, defining a suitable 
square probability matrixM, which is calculated for both spin orientations, up (f) and down 
(|), and for the two magnetic configurations, APC and PC. The component M a p of the 
matrix above, represent the probability that an electron which goes into piece a will come 
from piece j3. The latter takes into account conservation laws for scattering processes at the 
interface. 

Defining [A] as being a column matrix whose components are the vectors A( a ) (a = 
1, ...,n), and integrating (||) backwards in time, we can write the following matrix equation 
for A( a ): 

[A]=M{E[A] + [V]}, (5) 

where [V] (a column matrix) and E (a diagonal square matrix) are spin and configuration 
dependent. The components of [V] and E have the form 

V a = e~^ T v(t', a)e t 'l T di!, E af3 = e'^S^ . 
Jo 

Inserting the solution of ([5]) in (^), and then in ([I]), we obtain the conductivity tensor. 
Finally, by tensor inversion of the latter, we obtain the resistivity. The magnetoresistance 
(Ap/p) is defined here as 

Ap/p = On -ph)/(pu + ptt) > ( 6 ) 

where p^ and p^ are the resistivities of the sample for APC and PC, respectively. We note 
that the magnetoresistance considered here is due to the internal magnetic field, which for a 
given sample is a constant. To vary % = uj c t in relation implies measurements in several 
samples with different purities and internal fields (neglecting variations with temperature). 

We consider several topologies and connectivities induced by different values of rj. The 
limit cases of high (T — ► 1) and low (T — > 0) transmissions, are shown in Fig. ([[]). The 
examples include open orbits for majority (minority) spin and closed ones for minority 
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(majority), and closed orbits for the two spins. In all calculations we assumed a constant 
and isotropic effective mass of carriers; furthermore, we consider the mean collision time r, as 
being the same across the sample. The spin dependent behavior is then entirely due to Fermi 
surface effects. We vary the height and the width of the barrier at the spacer, assuming values 
given in the literature for some metallic multilayers. The magnetic potential, depending on 
the exchange splitting, is assumed to be of the order of eV, as in the case of ironi. 

Some symmetry considerations are in order. In the APC configuration, when going from 
one magnetic layer to the next, majority and minority spins are interchanged. Thus, elec- 
trons with spin f and j have the same density and topology in the limit cases depicted in 
Fig. [1]. Antiparallel magnetic fields in neighboring layers drive the electrons in orbits with 
different curvature, and one may get compensated-like orbits in some cases. This effect is 
clear when we consider closed orbits for the two spin states (see Fig. 0(c)), but it can happen 
in other cases in which open orbits break down in identical closed orbits, in the high trans- 
mission limit (the case shown in Fig. |T](b) is an example of this effect). Compensated-like 
orbits highly enhance the magnetoresistance in both directions x and y in a nearly isotropic 
fashion. Compensation between the two spin states for the APC case, yields a quadratic 
increase of the resistivity with the field. On the other hand, for the PC configuration, all 
closed orbits are electron-like, leading to saturation of the magnetoresistance. The net effect 
for Ap/p is a huge magnetoresistance effect, which is shown in Fig.|2] as a function of the 
dimensionless parameter x — UJ c T - 

We have also made calculations that include open orbits, with open orbit orientations 
perpendicular (Fig. [l|(a)) and parallel (Fig. |l](b)) to the interface. Figure |] shows some 
particular cases calculated. All the curves are plots of Ap/p as a function of x- The 
magnetoresistance is anisotropic when an open orbit is present. In same cases we may get 
Ap/p < values, i.e. inverse magnetoresistance effect^ (see left side of Fig.|3|). 

In summary, diamagnetic effects are important due to high internal fields in ferromagnets. 
We have shown that compensation-like effects in the APC configuration yield high values 
for the magnetoresistance, when the scattering time t a is equal or longer than the cyclotron 
period T c . In general, when r] ^> 1, electrons have enough time to probe the Fermi surface 
topology, before being scattered at the interface. 

On the other hand, under particular conditions that lead to open orbits, one may get 
the inverse magnetoresistance effect. The latter phenomenon depends critically on the 
isochronous assumption for the orbit pieces that are connected by scattering at the in- 
terface. For a general distribution of scattering times, electron orbits will trace random 
walk trajectories, with the corresponding loss of coherence. 

Scattering from the interface between magnetic layers produces, in general, a mixing of 
orbits as in 'breakdown' phenomena, suggesting a mechanism for the GMR based exclusively 
on Fermi surface effects. Our results indicate that the so called Kohler's rule is not valid, 
since there are two relevant parameters, x an d tj, in the theory. The authors are currently 
extending the present calculation to include a finite spin-flip relaxation time, i.e. a finite 
spin diffusion length.0 

The authors acknowledge partial financial support from Conselho de Desenvolvimento 
Cientifico e Tecnologico (CNPq, Brazil) through projects # 143009/95-8 and # 301221/77-4. 
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FIGURES 



FIG. 1. Connectivities and orbit topologies considered in our calculation, in the limits T — > 
and T — > 1, for the transmission coefficient at the interface. Cases (a) and (b) correspond to 
7) = 1/2, and (c) to rj = 1. We depict both magnetic arrangements, with the APC and PC 
configurations. In (a), orbits for majority and minority spins that are closed at T — > 0, breakdown 
to open ones in the direction perpendicular to the interface, when T — > 1. In (b), we have, at T — > 0, 
closed orbits for minority spin and open ones for majority. They breakdown, at T — > 1, to different 
topologies depending on the magnetic configuration: for the APC setup, we get compensated-like 
closed orbits; for the PC one, we get an open orbit for minority spin plus closed ones for majority. 
Case (c) is the same as (a) at T — > 0, but now the different value of r] (rj = 1) implies a different 
connectivity at T — > 1. In this latter we get exclusively closed orbits, which are compensated only 
for the APC case. The thick line for the PC setup means that orbits are retraced. 

FIG. 2. Magnetoresistance as a function of uj c t for case (c) of Fig.|], near the high transmission 
regime (T ~ 1). The parameters are indicated in the figure (mass of the carriers m, Fermi energy 
Ep, height of the interface barrier Vb, and width of the spacer d). The spin splitting V\ is also 
indicated in the inset of the upper figure. The magnetoresistance is nearly isotropic, as shown in 
the figure for the x and y directions. Note that the huge GMR effect is due to compensation in 
the APC setup, and is essentially independent of the splitting and the transmission coefficient at 
the barrier. 

FIG. 3. Magnetoresistance as a function of uj c t for the cases (a) and (b) of Fig.[l], at fixed 
exchange splitting and for several values of the interface barrier Vb, which are shown in the upper 
right figure. Other parameters are the same as in the previous figure. Both cases are anisotropic, 
due to the presence of open orbits. Case (a) displays the inverse magnetoresistance effect along 
the y direction, for values of the parameters near the case T ~ 1. 
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